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ABSTRACT 



Using a new XMM-Newton observation, we have characterized the X-ray prop- 
erties of the middle-aged radio-quiet 7-ray pulsar J0357-I-3205 (named Morla) and 
its tail. The X-ray emission from the pulsar is consistent with a magnetospheric 
non-thermal origin plus a thermal emission from a hot spot (or hot spots). The 
lack of a thermal component from the whole surface makes Morla the coldest 
neutron star in its age range. We found marginal evidence for a double-peaked 
modulation of the X-ray emission. The study of the 9'-long tail confirmed the 
lack of extended emission near the pulsar itself. The tail shows a very asymmet- 
ric brightness profile and its spectrum lacks any spatial variation. We found the 
nebular emission to be inconsistent with a classical bow-shock, ram-pressure dom- 
inated pulsar wind nebula. We propose thermal bremsstrahlung as an alternative 
mechanism for Morla's tail emission. In this scenario, the tail emission comes 
from the shocked interstellar medium (ISM) material heated up to X-ray temper- 
atures. This can fully explain the peculiar features of the tail, assuming a hot, 
moderately dense interstellar medium around the pulsar. For a bremsstrahlung- 
emitting tail, we can estimate the pulsar distance to be between 300 and 900 
pc. A pulsar velocity of ~1900 km is required - which would make Morla the 
pulsar with the largest velocity - and high inclination angles (>70°) are preferred. 
We propose Morla's nebula as the first example of a new "turtle's tail" class of 
thermally-emitting nebulae associated to high velocity pulsars. 

Subject headings: Stars: neutron — Pulsars: general — Pulsars: individual (PSR 
J0357-F3205) — X-rays: stars 
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Introduction 



The dramatic increase in the number of known gamma-ray pulsars since the launch 
of the Fermi Gamma-ray Space Telescope offers the first opportunity to study a sizeable 



population of t hese 



Atwood et al. 



fAbdo et al. 



2009) 



2010a 



ligh-energy sources. The Fermi Large Area Telescope (LAT, 
) has discovered pulsed 7-ray sign als from more than 100 objects 



LAT collaboration 



in preparation! ), revolutionizing our view of 



them and marking t he birth of new high-energy pul sar sub-families, suc h as millisecond 



flAbdo et al. 



2009al ) and radio-quiet 7-ray pulsars (lAbdo et al. 



2009bf ). that are as 



numerous as the classical one. The wealth of detections confirms the importance of the 
7-ray channel in the overall energy budget of rotation-powered pulsars and points to 



emission models in which the 7-ray 
open-field lines (outer gap/slot gap, 



production o ccurs i n the outer magnetosphere a. 



Cheng et al. 



(119861 ) 



Harding fc Muslimov 



ong 



toom and 



paves the way for our understanding of the three-dimensional structure and electrodynamics 
of neutron star magnetospheres. 

In a sense, the most important objects to constrain pulsar models are the "extreme" 
ones, accounting for the tails of the population distribution in energetics, age and magnetic 
field. In this respect, PSR J0357+3205 is one of the most interesting pulsars discovered 
by the LAT. J0357+3205, which we name "Morla" (cf. "The Neverending Story" , M. 
Ende) for its slow rotation among Fermi-LAT pulsars (P=444 ms) , was discovered in a 



blind search during the first three months of Fermi-LAT scanning ( lAbdo et al. 



2009cf ). It 



is located off the Galactic plane, at a latitude of b~-16°. While Morla is not the oldest 
rotation-powe red pulsar (r^=5.4 x 10^ yr) - with an age similar to the well-known "Three 



Musketeers" (jPe Luca et al. 



20051 ) - it is the slowest rotator among Fermi pulsars. Its 



spin-down luminosity, as low as Erot = 5.9 x 10^^ erg s ^, makes it one of the non-recycled 
7-ray pulsars with the smallest rotational energy loss detected so far. This suggests that 
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PSR J0357+ 3205 is quite nearby: by sca ling; its 7-ray flux using the "pseudo-distance" 



relation from 



Saz Parkinson et al. 



( I2OIOI ). we obtain a distance of ~500 pc, making Morla 



a natural target for X-ray observations. 



De Luca et al. 



mm found the X-ray 



Exploiting a joint Chandra and A/^O^IO program, 
counterpart of the 7-ray pulsar to be located at R.A.(J2000) = 03'*57™52^32, Dec.(J2000) 
= +32°05'20.6" (0.25" error radius, 68% confidence level). The pulsar X-ray emission 
was found to be consistent with a purely magnetospheric, non-thermal origin leading to 
a value of the 7- t o-X fl ux ratio of ~1000, well below the mean for radio-quiet pulsars 



-9000, 



Marelli 



( 120121 )). Unfortunately, the poor statistics precluded a clear detection of 



a thermal component in the neutron star spectrum as well as a full characterization of its 
column density and distance. Nevertheless, the Chandra data unveiled the presence of an 
extended structure protruding from the pulsar and extending >9' in length, which would 
correspond to ~1.3 pc for a pulsar distance of 500 pc, and ~1.5' in width. The "tail" is 
apparently detached from the pulsar, without detection of nebular counts until 50" from 
the pulsar; it shows an unusual, asymmetric morphology, with the brightnest part far from 
the pulsar. NO AO optical images, taken in V and K bands, show no hints of correlated. 



diffuse emission. No extended source was detect ed either in pub 



The lack of any Ha emission around the pulsar (jPe Luca et al. 
full ionized, likely hot (> 10^ K), ISM. 



ic radio or infra red surveys. 



in preparation 



) points to a 



Recently, using a new Chandra observation. 



De Luca et al. 



(lin preparation! ) detected the 



pulsar's proper motionof 165±30 mas yr ^, in the direction away from the tail and parallel 



to the Galactic plane. Assuming a pulsar distance of 500 pc, this results in a lower 



velocity of 390 km s ^, making it a relatively high- velocity pulsar ( iManchester et al. 



imit 



2OO5I). 



Many extended tails of X-ray emi ssion have been d i scover ed, associated with a number 



of rotation powered pulsars (e.g. 



Kargaltsev et al 



ram-pressure dominated pulsar wind nebula (see 



fl2008h) within t 



le framework of a 



Gaensler fc Slane 



(|2006[ ) for a review) 
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However, the unusual phenomenology of the feature we observed is barely consistent with 
such a picture. There is no hint of diffuse emission in the pulsar surroundings, where bright 
emission from the wind termination shock should be seen (as observed in other known 
cases). Moreover, the asymmetric brightness profile of the tail, with its sharp north-eastern 
edge and its broad maximum at large distance from the pulsar is remarkably different from 
any other observed structure (with the possible excepti on of the peculi a r feat ure associated 



with PSR B2224+65, powering the Guitar nebula - see 



Hui fc Becker 



(120071 )1 The mere 



existenc e of the tail is also prob 



(see e.g. 



Jager fc Djannati-Atai 



ematic for a pulsar with such a low Erot as J0357+3205 



In order to assess the nature of the tail and its relationship with the pulsar, as well 
as to better constrain the pulsar emission properties and physics, we obtained a deep 
observation of the pulsar and its peculiar tail with XMM-Newton. Exploiting its high 
spectral resolution, we were able to disentangle the components of the pulsar spectrum 
and constrain its distance based on the absorption column density. A complete study of 
the timing behavior of Morla in X-rays was carried out as well. Finally, we studied the 
properties of the tail, in order to characterize the mechanisms responsible for its emission. 



Observations and Data Reduction 



Our deep XMM-Newton observation o f Morla started on 20 11 September 15 at 02:37:18 



UT and lasted 111.3 ks. The PN camera ( Struder et al. 



200 ih of the EPIC instrument 



was operating in Large Window mode (tim e resolution of 47.7 ms over a 14' x 27' field 



of view (FOV)), while the MOS detectors (iTurner et al 



2001 



) were set in full frame 



mode (2.6 s time resolution on a 15' radius FOV). The thin optical filter was used for the 
PN camera while we chose to use the medium filter for the MOS detectors due to the 
presence of moderately bright (m/j ~9) stars. We used the XMM-Newton Science Analysis 
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Software (SAS) vll.O. After standard data processing, using the ep proc and emproc too. 



and sc reening for high particle background time intervals (following 



De Luca fc Molendi 



(I2OO4J )). the good, dead-time corrected exposure time was 98.5 ks for the PN and 108.3 ks 



for the two MOS detectors. In order to fully cha racterize both the pu lsar and the nebula 



we also used the three available Chandra/ ACIS ( iGarmire et al 



2OO3I ) observations of the 



field, performed on 2009 October 25, 26, and 2011 December 25, lasting 2 9.5, 47.1, and 



29.0 k s res pectively (these data sets were i ncluded in the investigations by 



torn and 



De Luca et al. 



De Luca et al. 



(lin preparation! )). We retrieved "level 2" data from the Chandra 
Science Archive and used the Chandra Interactive Analysis of Observation (CIAO) software 
v.4.2. We also combined the Chandra spectra, response matrices and effective area files 



using the mathpha, addarf and addrmf tools. 



2.1. Imaging and Source Detection 

Figure [1] shows the 0.3-6 keV exposure-corrected XMM-Newton FOV image, where the 
PN and two MOS images were added with the ximage tool. We applied a Gaussian filter 
with a kernel radius of 5". The image clearly shows the straylight effect, probably due to 
the bright source X Persei that dominates the X-ray sky ~70' to the southwest of Morla, 
with a flux of ~10~^ erg cm^^ s~^. X-ray straylight in EPIC is produced by photons which 
experience single reflection by the mirror hyperbolae and which reach the sensitive area 
of the camera and produce near-annular structures in the image. In order to study the 
background sources in the XMM-Newton observation, we evaluated the spectrum of our 
str aylight annuli. This is consistent, a t low energies, with the X Persei spectrum found 



by 



La Palombara fc Mereghetti 



(l2007l ). thus we added and froze these components in the 



spectral model of sources affected by this phenomenon. 

Source detection using maximum likelihood fitting was done simultaneously on each 
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of the EPIC-PN, MOSl, and M0S2 using the SAS tool edetect_chain. This tool runs 
on the event lists and invokes several other SAS tools to produce background, sensitivity, 
and vignetting-corrected exposure maps. A likelihood threshold of 10 was used for source 
detection, corresponding to a significance level of 3.6a. Figure [1] shows the brightest sources 
in the PN field of view; they were analyzed in order to produce X-ray spectra and find 
optical counterparts. The first 7 objects were also observed by Chandra so that we used 
both the telescope datasets to obtain radial profiles and spectra. 

First, we analyzed their radial brightness profile and compared them with the XMM-Newton 
and Chandra theoretical Point Spread Functions (PSFs). We found sources 8 and 14 to 
be the superimposition of two different point-like sources. Using the 2011 Chandra data, 
source 7 is incompatible with a simple point-like source. This is in agreement with results 
from XMM-Newton data - though heavily affected by straylight photons - of a point-like 
source surrounded by extended emission. The NASA/IPAC Extragalactic Database (NED 
j) reports one radio source compatible with the point-like emission, identifying it as an 
Active Galactic Nucleus (AGN) in a cluster. 

Each spectrum was then fitted either with an absorbed powerlaw, well-suited for AGNs 
as well as pulsars and an absorbed apec, well-suited for stellar coronae. Only source 9 
requires a double apec, quite common for stellar coronae; for this source we thus cannot 
exclude the superimposition of two different stars closer than ~2". Three sources (11, 13, 
and 14) cannot be described by the simple models considered here; deeper studies on the 
serendipitous sources are beyond the aim of this paper so they were excluded from our 
analysis. 

Optical counterparts were searched in the deep (4hr 26min exposure) optical and near- 
infrared images collected with the 4m Mayall Telescope on 2009 November 10, using the 



http:/ /ned. ipac.caltech.edu/ 
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large-field (36'x36') Mosaic CCD Imager (jjacoby et al. 



19981 ). All the results of our X-ray 



and optical analyses on background sources are reported in Table [TJ 
We found 2 sources (8 and 9) that are compatible with a star interpretation both in the 
X-ray and optical bands. These two objects are likely located closer to us than Morla, given 
their lower value of the absorption column. Deeper optical observations with different filters 
would be required to find their distances from the optical absorption. All the remaining 
objects are extragalactic, while source 17 presents a ve ry high value of the X -ray to optical 



20021 ). It is possible 



fiux ratio. This can be described as an obscured AGN (iMainieri et al. 
that the emission actually comes from two different sources (thus lowering the X-ray to 
optical flux ratio), though the low statistics prevent us from studying the brightness proflle 
of the source. 

All the AGNs were fltted together, linking their Nh to assess the average Galactic absorption 
column. The resulting Galactic valu e N^' = 2.10±0.09 x 10^^ cm^ ^ is in good agre e ment 



with the theoretical one predicted by 



Dickey fc Lockman 



fll990h and 



Kalberla et al. 



fl2005h 



for objects at the limit of the Galaxy (~500 pc). 



The Pulsar 



3.1. Pulsar Spectral Analysis 



In order to optimize the pulsar extraction radius we evaluated the brightness proflle 
of the source to maximize both the signal-to-noise ratio (S/N) and the number of counts. 
The PSF of the EPIC-PN camer a onboard X MM-Newton is best described by an off-axis. 



energy-dependent King function (IRead 



2004 ). The full width half maximum of the PSF 



for an on-axis source at 1.5 keV is typically less than 12.5" for the PN camera and 4.4" 
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Table 1: X-ray and optical analysis of background sources 



Object 








Coordinates 


specx 


(10^^ cm ^) 




type 


1 


03 


57 


47 


50, +32:03:45. 10 


pow 


'X o o 1 r\ r\r\ 

1.83+0.20 


5.0 


AGN 


2 


03 


57 


50 


95, +32:01:47. 70 


pow/apec 


9 1 1 +1.20 /I oc+0.87 


6.5 


AGN 


3 


03 


57 


53 


90, +32:00:58. 50 


pow/apec 


l-02^SiV0-66l°:^« 


9.1 


AGN 


4 


03 


58 


21 


95, +31:59:11. 65 


pow 


1.87+0.22 


1.8 


AGN 


5 


03 


58 


15 


25, +32:00:17.35 


pow 


-I n-i +n fi6 

-'-■^-'--0.58 


1.0 


AGN 


6 


03 


58 


00 


29, +31:58:01. 90 


pow/apec 




1.5 


AGN 


7 


03 


57 


39 


or" iOO/^i r\r\ 

35, +32:01:03. 00 


pow/apec 


c\ o 1 r\ r\ 1 A o/^ -4-0. oo 

2.68±0.60/1.20_o;34 


9.2 


A AT 

AGN 


8a 


03 


57 


12 


7,5+32:00:50.10 


apec 


r>r+0 76 
-■-•^^-0.68 


0.0056 


btar 


8b 


03 


57 


12 


95, +32:00:42. 65 


apec 


^ , rv-l-n 25 


0.0028 


Star 


9 


03 


57 


07 


35, +32:02:59. 65 


2apec 


o.2i+°:}2 


0.0060 


Star 


10 


03 


57 


46 


30,+31:57:24.50 


pow 


1.98t°J^ 


3.6 


AGN 


12 


03 


57 


20 


20,+32:03:14.95 


pow/apec 


1.69+°|^/0.98+°:^^ 


0.12 


AGN 


15 


03 


58 


40 


80,+32:02:27.35 


pow 


1 O7+0.46 
-0.50 


2.8 


AGN 


16 


03 


58 


20 


35,+32:08:24.85 


pow 


i.83t}:?I 


1.0 


AGN 


17 


03 


58 


17 


45,+32:08:50.20 


pow/apec 


1.73+0.72/0.941°:^? 


400 


AGN 



''We report the R.A., Dec. (J2000) coordinates of the optical counterpart, with a typical error of 0.3". If 
there was no optical counterpart, X-ray coordinates are shown, with a typical systematic error of 2" and Id 
statistical error of 



Note. — X-ray and optical parameters of the brightest sources in the PN field of view; the colors are 
in a logarithmic scale. All the errors are at a 90% confidence level. We use both the optical to X-ray fiux 
ratio and the X-ray spectral type to separate AGNs from stellar coronal emissions. PSF studies show that 
sources 8 and 14 result from the superimposition of two sources while source 9 can be a single source. The 
spectra of sources 11, 13, and 14 cann ot be described by a p owerlaw nor apec components. Source 17 could 



be associated with an obscured AGN (Mainieri et al. 



20021). 
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Fig. 1. — Combined 0.3-6 keV exposure-corrected XMM-Newton field of view, where tlie 
PN and two MOS images were added witli ximage tool. We applied a gaussian filter with 
a kernel radius of 5". The yellow crosses mark the serendipitous sources we analyzed to 
constrain the value of the Galactic absorption column (see Table [T]). In the lower-right part 
of the figure the straylight effect is apparent. 



De Luca et al. 



a 



for the two MPS dete ctorss q Our results are in agreement with those pubhshed by 



( 120111 ). based on Chandra data, with no excesses within ~ 30" of the 
pulsar, with the exception of its bright tail. We found the best extraction radius to be 18" , 
containing 2022, 515 and 540 counts in the 0.3-10 keV energy range in the PN and the two 
MOS detectors respectively, with a background contribution of ~10%. We generated an ad 
hoc response matrix and effective area files using the SAS tools rmf gen and arfgen. We 
also analyzed the Chandra spectra of the source using a 1.5" extraction radius (698 pulsar 
counts with a background contribution of 0.4%); we added the three different spectra with 
the tool mathpha. Response matrix and effective area files were generated using the CIAO 
script psextract. Due to the presence of the pulsar tail, backgrounds were evaluated from 
a circle on the same CCD as the pulsar and without serendipitous X-ray sources, with a 
radius of 10" for Chandra and 50" for the XMM telescopes. 

For the pulsar emission, we tried a simple powerlaw (pow), the combination of a powerlaw 
and a blackbody (bb) as well as the combination of a powerlaw and a magnetized neutron 
star atmosphere model (nsa in XSpec - assuming a neutron star with a radius of 13 km, 
mass of 1.4 Mq and a surface magnetic field of 10^^ G). While the simple powerlaw gives a 
fit only marginally acceptable (reduced chi square x^=l-22, 134 degrees of freedom (dof)), 
both the composite models yield better fits (nsa-|-pow: x^=0.96, 131 dof; bb+pow x^=0.97, 
132 dof). An F-test performed comparing a simple powerlaw with the composite spectra 
results in a chance probability of less than 10^^, pointing to a significant improvement 
from adding a thermal component. The pulsar emission is well described by a powerlaw 
component with a photon index of F = 2.28]l^oJg (90% confidence errors) and a blackbody 
component with a temperature of kT = 94.0^9^2^ eV and an emitting radius at R = 
0.45^Q;'^g(i50o km (where d^oo is the distance of the pulsar in units of 500 pc), absorbed 



^http: / /xmm.esac.esa.int / external / xmm_user_support / documentation / sas_usg/USG/ 
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by a column density of Nh = 1.20lo'45 x 10^^ cm~^. The 0.3-10 keV non-thermal pulsar 
luminosity is Lnth = ^■^'^-o^l^loo ^ 10^° ^rg s~^, while the thermal one Lnth = l-14lo]2o^5oo 
X 10'^'^ erg s~^. As expected, the nsa model yields a lower temperature of 41.ll;g2 eV and 
a larger emitting radius of R = S.TTl'^logrfsoo km, and a powerlaw component with T = 
2.21±0.19, absorbed by a column density of Nh = lAdto^l ^ 10^^ cm~^ Both results ar e 



i n agr eement with the upper limits on the thermal components found by 



De Luca et al. 



( 1201 ll ). Figure [2] shows the results of the fit for the two different spectral models. 



3.2. Pulsar Timing Analysis 

In order to carry out a proper timing analysis with EPIC data, a correct photon arrival 
time must be assigned to each event detected on board. The transformation from readout 
sequences by the EPIC camera to photon arrival times of each photon is performed by the 
EPEA (European Photon Event Analyser, 



Kuster et al. 



(119991)). A problem in the EPEA 



can produce time jumps in some observations, which have to be corrected. In order to 
identify any possible remaining time jumps, one must neglect both the barycentering and 
the randomization of photon arrival times. The time At between successive events is then 
calculated and divided by the Frame Time (FT) of the relevant mode (47.66384 ms for 
PN Large Window mode). A time jump is shown to exist when At/FT is different from 
an integer by a quantity larger than a tolerance parameter. Only those time jumps which 
hap pen to be an in teger multiple of the relevant FT would be missed by this method (see 



e.g. 



Marelli 



mm ). 

Our XMM-Newton data show that a time jump happened in all the three instruments 
between 2011-09-15 16:04:45.822 and 16:10:14.310 UTC, at about the half of the observation. 
Such event makes it more difficult to find the pulsation of Morla, creating a possible 
misalignment in the pulsar phase assigned to the photons before and after the jump. For 
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Fig. 2. — PN and MOS spectra of the Morla pulsar. Pattern 0-4 PN events and pattern 
0-12 MOS events have been selected among photons within 18" of the target position. The 
spectra are rebinned in order to have at least 25 counts per bin. The black, red, and green 
curves show the data from the PN and two MOS detectors, respectively. Also shown are 
the spectral fits: magenta curves show the thermal components, while orange curves the 
non-thermal ones. The data in the left panel are fitted with an absorbed powerlaw plus 
a blackbody while the data in the right panel use an absorbed powerlaw plus neutron star 
atmosphere (nsa) model. The lower panels show the residuals, in units of standard deviations, 
with 1 standard deviation error bars. 
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this reason, we performed careful timing studies in the two subsets, as well as in the total 
data set. 

We evaluated the optimal selection of pulsar events by maximizing the signal-to-noise 
ratio in the 0.2-6 keV energy range as a function of the source extraction region. The 
optimal choice turned out to be a source extraction radius of 18" for the PN camera. We 
selected only pattern photons in the range 0.2-0.35 keV in order to reduce the background 
contribution due to low-energy electronic noise, while in the 0.35-6 keV energy range, 
we used pattern 0-4 photons. Photon arrival times were converted to the solar system 
barycenter using the SAS tool bary cen. According to the accurate Fermi-LAT ephemeris 



( ILAT collaboration 



in preparation! ) . a period of P = 0.44410548661(9) s is expected at the 
epoch of the XMM observation (the uncertainty on the last digit is computed by propagating 
the 3(7 error of the Fermi-LAT timing solution). The uncertainty of ~ 10"^° s is much 
smaller than the intrinsic resolution of the XMM-Newton dataset (P^/2At ~ 8.4 x 10~^ 
s). Thus, if the spin-down of the pulsar is stable, a single trial folding is appropriate. 
Folding using the full LAT ephemeris yields marginal evidence for a modulation of the 
X-ray emission (see Figure [3]). A Pearson's test on the folded light curve (10 phase bins) 
in the 0.2-6 keV energy range results in a probability of 5 x 10~^ of sampling a uniform 
distribution. Due to the low signal modulation, a study of the phase-misalignment before 
and after the time jump was not possible. To evaluate the significance of the detection 
independently of the ph ase binning chosen, w e computed the corresponding value 



using the Rayleigh test ( IBuccheri et al. 



19831 ) for 1 and 2 harmonics. This exercise yields 



probabilities as small as 5 x 10^^ and 3 x 10^^, respectively, for a chance fluctuation. Such 
a test confirms the hint of a modulation at the pulsar frequency (and its first harmonic) in 
the XMM-Newton data, with a misalignment between the X-ray and gamma-ray peaks (see 
Figure [3l). However, the low number of counts prevented us from disentangling the thermal 
from the non-thermal pulsation or from performing a phase-resolved spectroscopy analysis. 
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Fig. 3. — EPIC/PN folded light curves in different energy ranges using photons within a 
18" radius of the Chandra position. X-ray photon phases were computed by folding at the 
expected period, based on the most recent available LAT ephemeris: the pulsar period at 
the start of the observation is P = 0.44410548661(9)s and the P contribution is taken into 
account. We note that the XMM-Newton time slip (discussed in section 13. 2p could affect our 
results. Pattern events have been selected in the 0.2-0.35 keV energy range while Pattern 
< 4 have been used in the 0.35-6 keV range. For each energy band, the relative contributions 
of all the spectral components are showen. On average the background accounts for ~15% 
of the counts. 80% of the 0.8-6 keV counts come from the non-thermal pulsar component 



while 50% of the 0.2-0.8 keV counts come 
LAT light curve of the Morla pulsar from 
light curves have been aligned in phase. 



rom the thermal on e. The lower panel shows the 



AbdoetaL 



fl2010al ) to which the XMM-Newton 
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4. The Tail 

4.1. Spatial Analysis 

The 2009 Chandra observations revealed an extended, bright nebula protruding from 
the pulsar, extending ~9' in length. Unfortunately, the 2011 Chandra observation does not 
cover the entire nebular emission so it cannot be used. The extended feature can be easily 
studied in our XMM-Newton observation due to its lower particle background and higher 
spectral resolution, compared to the Chandra one. On the other hand, the worse spatial 
resolution of XMM-Newton prevented us from significantly constraining the properties of 
the nebula near the pulsar. We therefore concentrate on a higher-scale analysis of the 
extended feature. 

First, we produced brightness profiles of the nebula along the pulsar motion (APM) 
and orthogonal to it (0PM). In order to maximize the statistical significance of each bin, 
we carefully evaluated the bin size and the width of extraction boxes, choosing a scale of 
90"x50" APM and 450" x 12. 5" 0PM (see Figure i]). For each instrument, we excluded 
the emission region surrounding the pulsar, as well as all the serendipitous sources and 
CCD gaps, according to our source PSF and calibration studies. Area calculations for each 
bin take into account the excluded regions. The presence of straylight annuli, PN CCD 
gaps, and numerous serendipitous sources in the faintest part of the nebula, toward the 
south-west, forced us to choose ad hoc boxes in that region, in order to exclude all the 
contaminants. 

The brightness profile for both XMM-Newton and Chandra are shown in figure O In order 
to take into account the different angular resolution of the XMM-Newton instruments, 
the profiles were convolved with a Lorentzian kernel of F = 7.125", corresponding to a 
full-width half maximum (FWHM) of the distribution of 6.6", the nominal FWHM of the 
XMM-Newton mirrors. Figure [5] shows the smoothed Chandra and XMM-Newton brightness 
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profiles. While minor differences are found in the count rate profiles shown in the figure, 
taking into account both the different (energy-dependent) effective areas, the flux profiles 
are compatible. The XMM-Newton data show a hint of extended emission, at brightness 
below the Chandra detection limit, near the pulsar itself. Such a difference is possibly due 
to the higher XMM-Newton sensitivity, that highlights the fainter part of the nebula, right 
behind the pulsar. No variations in the tail profile are detected in the available observations. 
Along the pulsar motion, the tail emerges from the background within 50" of the pulsar 
and its flux increases gradually, reaching a flat maximum between 4' and 7'. Its decrease 
is faster, though the profile is affected by the two serendipitous AGN (see section 12. ip and 
it fades below the background level at about 9' from Morla. In the direction orthogonal 
to the main axis (0PM), the tail profile is clearly asymmetric: it shows a fast increase in 
the northeast direction (rising to the maximum within 30"), reaches a sharp maximum, in 
broad agreement with the neutron star proper motion direction, then the decrease is slower, 
fading at about 1-1.5' in the nearest part of the nebula and 2' in the furthest. Such a 
behavior recalls the projection of a cone, with the apex corresponding to the pulsar position 
(or somewhere in the vicinity of the pulsar). 

Brightness profiles of the tail in the soft (0.3-1.5 keV) and hard (1.5-6 keV) energy bands 
show no significant discrepancies, with the possible exception of the region immediately 
next to the pulsar (<50"), where a hint of a softening in the spectrum is seen. 

4.2. Spectral Analysis 

While in the previous Chandra observations the spectral analysis of the tail emission 
was hampered by the low signal-to-noise ratio (S/N) - ~57% of the total counts coming 
from the background, - the exceptional sensitivity of XMM-Newton allows for a deep study 
of the nebular emission, with a background contribution of ~35% . 
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Fig. 4. — 0.3-6 keV exposure-corrected XMM-Newton view of the tail, where the PN and 
two MOS images were added with the ximage tool. Left panel: The green boxes mark 
regions used for the brightness profile study along the pulsar proper motion (indicated with 
an arrow, labeled APM) while white ones mark regions used for the brightness profile study 
orthogonal to the pulsar proper motion (indicated with an arrow, labeled 0PM). Yellow 
boxes mark the ad-hoc regions used in the region contaminated by straylight annuli. The 
barred green circles mark regions excluded from the analyses due to the presence of point- 
like sources. Right panel : The red polygon mark the previously defined regions in which 
the excess due to the presence of the nebula is significant at 5a or more. Following the red 
polygon, we traced in magenta a plausible shape of the tail while the cyan ellipse is the 
extraction region used for our spectral analyses on the tail. The long green arrow indicates 
the direction of the pulsar proper motion, while the two shorter arrows represent the 1 sigma 
confidence errors on this direction. 
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Fig. 5. — Background-subtracted surface brightness profiles of tlie tail along its main axis 
(APM, left panel) and orthogonal to it (0PM, right panel). The FN plus MOS profiles in 
the 0.3-6, 0.3-1.5 and 1.5-6 keV energy ranges are marked, respectively, with black, red, and 
green. The cyan line shows the smoothed 0.3-6 keV C/ian(ira brightness profile. This last has 
been convolved with a Lorentzian kernel of F = 7.125" , corresponding to a Full width at half 
maximum (FWHM) of the distribution of 6.6", the nominal FWHM of the XMM-Newton 
mirrors. The integral of the Chandra profile has also been normalized to the integral of the 
0.3-6 keV XMM-Newton one in order to show any possible discrepancy. The blue vertical 
lines mark the position of the pulsar. Sources 4 and 5, indicated with an arrow, have been 
removed, possibly enhancing the decrease in brightness. The fiat maximum along the main 
axis is clearly visible between ~250" and 500" as well as the asymmetrical shape of the tail 
in the direction orthogonal to it. The profiles in the soft and hard X-ray bands are almost 
the same, with a possible exception right behind the pulsar, where the soft band dominates. 
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In order to maximize the S/N, we extracted the total spectrum from the brightest part of 
the nebula, an elliptical region centered at R.A., Dec. (J2000) 03:58:06.59, +32:02:05.42, a 
semi-major axis of 4' and a semi-minor axis of 50". We excluded from the analysis region 
a circle with 20" radius around source 5 and the gaps in the CCDs. We extracted the 
background from source-free circular regions uncontaminated by the extended emission 
and located on the same CCD as the tail. We generated the response and effective area 
files using the SAS tools r mf gen and arfgen. Chandra data were treated in the same 



way as m 



De Luca et al. 



( 120111 ) and the resulting spectrum was fitted together with the 
XMM-Newton ones. The spectrum of the nebula is well described {xl = 1.06, 162 dof) by a 
non-thermal emission model (powerlaw with a photon index F = 2.07±0.08), absorbed by 
a column Nh = (2.61±0.23) x 10^^ cm~^. Figure [6] shows confidence contours for Nh and 
F for the diffuse feature and for the pulsar. The value of obtained for the tail seems 
to be in conflict (at the ~ 3a level) with both the one obtained for the pulsar, as well as 
with the Galactic value, obtained by fitting the spectra of the serendipitous AGNs in the 
field (see section 12. ip . Such a discrepancy between the values of led us to consider 
alternative models. The spectrum of the tail is well described {xl = 1.06, 162 dof) by a 
thermal bremsstrahlung model with a temperature kT = 3.75]^q 4q keV, absorbed by the 
more realistic column Nh = (1.72+0.16) x 1 0^^ cm~^. Other similar thermal models such 



as mekal (hot diffuse gas 



Mewe et al. 



constant temperature iBorkowski et al. 



constant temperature 



Borkowski et al. 



P985 



), nei (coUisional plasma, non-equilibrium. 



(2001)), and pshock (plane-parallel shocked plasma. 



(1200 ll )) with redshift fit equally well the spectrum 



of the tail and require a very low metal abundance. Table [2] reports the results of the fit for 
all the models. 

In order to study the spatial-spectral variability of the nebula, we divided the elliptical 
extraction region using the two axes, then we fitted both the absorbed powerlaw and 
bremsstrahlung models (see Table [2]) for each semi-elliptical region. In this way we obtained 
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spectra from the semi-elliptical regions nearest to and furthest from the pulsar, and those 
at the south-west and north-east. No significant (3a) variation was detected down to ~0.2 
in r or ~1 keV in kT. 



5. Discussion 



5.1. The Pulsar 



The X-ray spectrum of the Morla pulsar is a combination of thermal radiation from 
hot spot(s) superimposed on a non-thermal power- law component. A straight estimate 
of t he NS polar cap size, b ased on a simple centered dipole magnetic field geometry (see 



e.g. 



De Luca et al. 



(120051 )) and a standard NS radius of 10 km, predicts in this case a 
polar cap radius of about 85 m. The observed emitting radius is about 450 m, significantly 
smaller than the entire surface of any reasonable NS but definitely bigger than the expected 
polar cap. This is true both using the blackbody model and the NS atmosphere model, 
although the latter yields a larger emitting area and a lower temperature. Such a large 



size of the spots has also been 



bund in other well-known pulsars, like PSR B0656+14 and 



PSR B1055-52 (iDe Luca et al. II2005I ). and it is still poorly understood. The luminosity 
associated with the Morla polar cap is 1.1 x 10'^° erg s~^, a factor of a few lower than the 
luminosity expected from polar caps heated by currents. However, it is consistent 
with the expectations from polar cap heating due to the bombardment by parti cles created 



only by inverse Compton scattered photons (iHarding &: Muslimov 



2001 



2002). PSR 



J0357+3205 is close to the death line for production of e"*"/ pairs by curvature radiation 
photons and this could explain the reduced polar cap heating. 



Kaspi et al. 



(|2006|)) 



While the majority of X-ray emitting isolated neutron stars (e.g. 
exhibit thermal emission coming from the whole NS surface, no such a component is 
required to fit the XMM-Newton spectrum. Using the blackbody model, for a NS radius 
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Fig. 6. — Confidence contours from spectroscopy on tlie pulsar and tlie tail. Black, red, and 
green lines respectively show the error ellipses at 68%, 90% and 99% confidence levels for the 
absorbing column and the photon index. The Galactic from the study of serendipitous 
sources (see section 2.1) is also shown with errors at 1, 2, and 3 a indicated by different 
shades of blue. It is apparent that the Galactic value of Nh and the value obtained by 
fitting the powerlaw spectrum on the tail are only marginally compatible. 
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Table 2: Parameters of different X-ray spectral models of the tail 



Model 




r 


kT Abund'^ 


Flux'' 








keV 




tot-pow 


1.06 


2.61±0.23 2.07±0.08 




6.65±0.40 



tot-brem 1.06 1.72±0.16 

tot-mek 1.13 1.73±0.16 

near-brem 1.15 1.76±0.14 
far-brem 

sw-brem 1.33 1.72±0.14 
ne-brem 



3.75 



+0.48 
-0.40 



q 7C+0.49 
'"-0.38 

A 79+0.89 
^- ' ^-0.69 

o oq+0.38 
'-'■'-'^-0.33 



<0.099 



4.01 



+0.64 
-0.52 



A OQ+0.55 
^•^•J-0.46 



4.95±0.30 
4.57±0.28 

9 90+0.20 
^•^'J-0.34 

2.05±0.16 
1.83±0.16 
2.48±0.18 



"Me tal abundances (He f i xed a t cosmic). The elements included are C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, Ni, 



Anders fc Grevesse 



mm 



see 

^Unabsorbed 0.3-10 keV flux. 

Note. — Here we report the best fits of models that can describe the nebular emission. The elliptical 
extraction region was chosen after an accurate study based on the signal-to-noise ratio. We used an absorbed 
powerlaw, bremsstrahlung and mekal, respectively. Thermal models obtain similar temperatures, absorption 
columns and abundances. The last four spectra describe the emission from the nearest and the furthest 
semi-elliptical nebular regions from the pulsar, the south-west and north-east ones, respectively, using a 
bremsstrahlung model to fit data. To describe the total tail emission, we simultaneously fitted the XMM- 
Newton spectra and the Chandra ones. For the different parts of the tail we used only XMM-Newton spectra 
(due to the different PSFs of the two telescopes) and wc simultaneously fitted the furthest with the nearest 
nebular spectra and the south-east with the north-west ones, chaining only the Nh values. All the errors 
are at a 90% confidence level. 
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of 10 km and a distance of 500 pc, the 3o" upper limit on the temperature is 4.4 x 10^ K. 
This constraint makes Morla the coldest NS in its age range, shghtly cooler than the coeval 



radio-quiet pulsar RX J1856-3754 (e.g. 



Sartore et al. 



(|2012[ )). Such a difference can be 



attributed t o differ ent magnetic fields, and thus different magneto-thermal evolution (e.g. 



Pons et al. 



(l2009f l). The non-thermal X-ray luminosity of Morla is in broad agreement 



with the relation betw een the X-ray 



spin-down lum inosity ( iPossenti et al. 



2002 



Kargaltsev et al. 



2008 



pulsars and t 


leir 


Marelli et al. 


2011 



Marelli 



I2OI2I). 



5.2. The Tail 



The morphology of the tail, its spatial extension - ~1.3 pc at 500 pc assuming 
no inclination with respect to the plane of the sky (z), - and the lack of any other 
related Galactic source led De Luca et al. 2011 to associate unambiguously the nebular 



emiss ion with the pulsar. 



(e.g. 



Kargaltsev fc Pavlov 



l ongate d X-ray tails coupled with pulsars are quite common 



( l2008al )) and are usually interpreted within the framewor k 



Gaensler fc Slane 



of bow -shock, ram-pressure-dominated pulsar wind nebulae (e.g. 
( I2OO6I )). When a pulsar moves supersonically, the shocked pulsar wind is expected to 
flow in a region downstream of the termination shock (the cavity in the interstellar 
medium created by the moving Neutron Star and its wind), confined by ram pressure. 
X-rays are produced by synchrotron emission from the wind particles accelerated at the 
termination shock, which is t ypically seen (if it can be reso lved) as the brightest part of the 



extended structure (see, e.g.. 



Kargaltsev fc Pavlov 



from magnetohydrodyn ar nic simulations 



Bucciantini et al. 



20051). 



De Luca et al. 



(iBucciantini 



(l2008br)\ . according to the 



2002 



expect ations 



Van der Swaluw 



2003 



in preparation! ) have detected and measured a 



very large proper motion for the pulsar: 165±30 mas yr ^ - 390 km s ^ across the plane 
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of the sky, at 500 pc - in the direction opposite to the tail, thus compatible with such a 
scenario. 

The tail shows no other similarities with other, synchrotron-e mitting nebulae (for a 



complete discussion about Morla's synchrotron-emitting model see 



De Luca et al. 



mm)- 



T he first problem arise s frorn energetic requirements for the emitting particles. According 



to 



Goldreich fc Julian 



(Il969[ ). the maximum potential drop between the pole and the light 
cylinder (in an aligned pulsar) is A$ = (3 i^^ot /2c) In the Morla case, the electrons 
are accelerated up to a maximum Lorentz factor 'jmax ~10^, which can be considered as 
an upper limit for the electrons injected in the pulsar wind nebula. The low Erot of Morla 
requ ires an ambient ma g netic field as high as ~50 fiG - while the mean Galactic value is ~1 



/iG (jjansson fc Farrar 



( I2OI2I ) and references therein). However, one must also take into 



account the magnetic field carried by the particles themselves. An independent estimate 



of the mag netic field in t 



le tai l can be obtained from the measured synchrotron surface 



brightness ( iPavlov et al. 



20031). Taking into account the uncertainty on the parameters 
we can conclude that a reasonable value of the magnetic field in the tail is in the range 
20 — 100 /iG, consistent with the value of ~50 fiG obtained before. 

The second problem is the lack of diffuse emission surrounding the pulsar, where 
the emission from the w'ind termination shock should be brightes t , as o b served in all the 



other known cases (e.g. 



Gaensler et al. 



McGowan et al. 



(120061); 



Kargaltsev et al. 



( 12004] ): 

( l2008l )). Such a shock could be unresolved by Chandra only in the case of an extremely 
large ambient density (several hundred atoms per cm^) and/or a pulsar speed higher than 
some thousands km s~^. Taking into account the measure of the proper motion, the 
distance of the pulsar should be larger than 1.3 kpc or the inclination of the motion should 
be higher than 70°. Such conditions imply a tail at least 3.5 pc long, so that it becomes 
even more difficult to account for the energy of the particles. This picture implies also that 
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a significant fraction of tlie point-like flux comes from a non-thermal termination shock, 
reducing dramatically the power-law component of the pulsar. While it is acceptable for 
the non-ther mal component of the pulsar to be a factor of some smaller than the value 



we obtained ( Marelli et al. 



2011 



Marelli 



20121 ). the lack of nebular emission around the 



pulsar, where the synchrotron emission should be stronger, remains unexplained. 

The brightness profile along the main axis of the tail is consistent with no emission 
within 30-50" of the pulsar, then it slowly increases to reach a broad peak at ~4'. This 
behavior is remarkably different from that observed in all the k nown synchrotron nebula e 



Kargaltsev et al. 



(120081)) 



that show peaks close to the position of their parent pulsar (e.g. 
where the acceleration of the wind particles is the most important. It is very hard to explain 
in this scenario why the emission peaks at ~1 pc from the pulsar. 

The asymmetric brightness profile of the tail in the direction perpendicular to its main 
axis is characterized by a sharp northeastern edge and a slow decay toward the south-west 
direction. The synchrotron emission from the nebula is only marginally dependent on the 
ISM density, so that ad-hoc and large variations in the ISM are demanded, in order to 
explain this strange profile. Even if possible, such a case is hardly believable. 
Furthermore, the inferred synchrotron cooling time implies a significant spatial-spectral 
evolutio n with respect to the dis tance from the pulsar, as observed in all synchrotron 



nebulae (IKargaltsev et al. 



20081 ) . Thanks to the unique throughput of XMM-Newton, we 



could perform a deep spatial-spectral study of the tail, analyzing separately the furthest 
and the nearest parts of the nebula. No significant differences in the spectral index are 
detected (see section 14. 2p . The lack of variation cannot be easily explained in terms of a 
synchrotron nebula. The lack of such a softening would require some particle re-acceleration 
mechanism within the tail itself. 

The excellent throughput of the XMM-Newton instruments has revealed another problem 
of the synchrotron nebula scenario, connected to the value of the column density {Nh ) 
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obtained by fitting the non-thermal model of the nebula. While it is marginally compatible 
with the value obtained for the pulsar, the inferred N^^ appears to be larger than the 
Galactic one in that direction, measured by fitting serendipitous AGNs in the field (see 
Figure [6]) . 

As reported in section \4.2\ a thermal bremsstrahlung model fits equally well the 
spectrum of the tail and produces an estimate of Nh in agreement with that of the pulsar 
and the Galactic one. In the bremsstrahlung scenario, the tail emission arises from the 
shocked ISM material heated up to X-ray temperatures. Under the assumption of a strong 
shock- wave and from the junction conditions at the shock front, we can estimate the speed 
of the pulsar directly from the temperature of the shocked ISM material at the head of the 



bow shock: Vpsr = \/VokT j {3fimp), where fi is the molecular weight of the ISM gas and 
we assumed an adiabatic index jad = 5/3. The temperature obtained by fitting the tail 
spectrum with a thermal bremsstrahlung model, kT ^ 3.75 keV, implies a pulsar speed 
Vpsr — 1376/i~^/^ km s~^. From the values of the abundances resulting from the spectral fit, 
we can assume that ^ ~ 0.54, so the speed of the pulsar is about 1900 km s^^, which would 
make it the largest ever observed in a pulsar. If we combine this result with the pulsar 
proper motion we can estimate the inclination angle i of the pulsar motion (and thus of its 
tail), i ~ acos(0.2 x rfsoo). Moreover, through the value of the photon flux measured using a 
thermal bremsstrahlung model, it is possible to estimate the number density of the shocked 
ISM gas as n ^ 0.86^50!/^ atoms cm~^. 

Thermal bremsstrahlung emission acts as a cooling mechanism. The time scale for 
cooling can be estimated as tcooi — {E thermal / J) 1 where Ethermai — nkT is the thermal 
energy per unit volume and J is the volume emissivity. This approximation gives 



tcooi — 6 X lO^T^/^^ngs) ^ yr, where qb is the frequency-averaged Gaunt 



typically ranges between 1.1 and 1.5, depending on the plasma temperature (jPadmanabhan 



actor that 



- 28 - 



2000|). From our estimates of the temperature and number density of the shocked ISM gas 
we expect tcooi — 10^ yr. Considering a pulsar proper motion of 165±30 mas yr^^ and an 
extension of the nebula of about 9', then the age of the tail (~3300 yr) is much lower than 
its cooling time. This result can explain why we have not measured a significant evolution 
of the temperature of the tail as a function of the distance from the pulsar. On the other 
hand, the shorter than expected length of the tail can be explained assuming an ad-hoc 
decrease of a ~3 factor in the number density n beyond the end of the tail - the volume 
emissivity of the thermal bremsstrahlung emission is J cc ra^T^/^. Such a strong dependence 
on the ISM density can explain the asymmetrical shape of the nebula in the direction 
orthogonal to its main axis: this could be a consequence of an inhomogeneity of the number 
density in the nebula, quite typical in the ISM. 

From the conservation of mass, energy, and momentum on the two sides of the shock 
front and from the values measured in the tail (e.g. the pre-shock temperature Tq and 
number density no), we can calculate the same physical quantities outside the nebula: 

n _ (7 + l)M^ 

no (7 + l) + (7-l)(M2-l) ^> 

T ^ [(7 + 1) + 27(M^ - 1)] [(7 + 1) + (7 - 1)(M^ - 1)] 

To (7 + 1)2M2 ^ ' 

where M = 1/ sin(a) is the Mach number and depends on the cone angle of the tail (a) 
and an adiabatic index jad = 5/3 is assumed for the ISM gas. In the strong shock- wave 
limit (M >> 1) no = n/4 ~ 0.22dr^oo'^ atoms cm"'^ while To depends on the Mach number 
and thus on the inclination angle i. From the geometry of the cone angle of the nebula, a 
plausible value of Tq is in the range of 10^-10^ K, for i > 60°. The pre-shock temperature is 



consi stent with that of the hot phase of t 



[e.g. 



Bland-Hawthorn fc Reynolds 



le ISM, which fills a large fraction of the Galaxy 



the requirement of a full ionized ISM by 



(I2OOO)). Such a scenario is also in ful l agreement with 



De Luca et al. 



(lin preparation! ). We note that 



the density required by our model is slightly higher than expected for the hot ISM: this 
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may result from the concerted action of massive stellar winds and supernova explosions. 



Yao fc Wang 



(I2OO5I ) report the presence of dense (10~^-1 atoms cm~^), hot (~10® K) 
envelopes of gas in our Galaxy. On the other hand, thermal models such as mekal, nei, 
and pshock point to a very low metallicity of the ISM around the pulsar, making more 
difficult to explain this dense envelope as the result of a supernova explosion or stellar 
winds. Future multi-wavelength observations of the field and more detailed Galactic gas 
models could better constrain and explain the ISM composition around Morla. 

A shock-wave scenario can easily explain the lack of diffuse emission surrounding the 
pulsar. Most of the energy in the pre-shock flow is carried by the ions - the kinetic energy 
in the streaming of the electrons is less than 1/2000 of the total - but the electrons are 
generally responsible for the cooling of the plasma: this means that the electron temperature 
(Te) is responsible for the dynamics of the shock and its emitted spectrum. The electrons 
must be heated by the ions before the emission becomes detectable and this process takes 
some time. Coulomb heating of the electrons behind the shock proceeds at a rate 

dft tgq 

where Tj is the ion temperature and teq = 1 .l{Te^'^ /n) s ~70 kyr is the equipartition time 
for a fully ionized plasma with a Coulomb logarithm of 30 (where Tg is in units of K and 
n in cm~^). If all the post-shock energy is in the ions (Te ^ Tj), then the increase in Tg 
due to coUisional heating follows: Te(co//) = 1.7 x {not^Vpg^.^Y'^^ K, where ^3 is the age 
of the tail in units of 10'^ yr, no is the pre-shock number density and Vpsrs is the pulsar 
speed in units of 10^ cm s~^. Thus Te{coll) < 6.6 x lO'^d^^Q^ K, not significantly different 
from the value observed in the tail. Moreover, Te{coll) < 5.5 x 10^ K within ~ 30" of the 
pulsar (corresponding to ~200-300 years, from the measure of the proper motion). This 
fully explains the lack of diffuse emission surrounding the pulsar. 

So far, no evidence of other pulsars characterized by a thermal bremsstrahlung emitting 
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tail has been reported. This could be explained by the peculiar conditions it requires, e.g. 
very fast pulsars and an hot and dense ISM. With its estimated velocity of 1900 km s~^, 



J0357+3205 is the fastest moving pu 



Pulsar Catalogue (IManchester et al. 



sar kn own. Among all the pulsars listed in the ATNF 



20051 ). only 2 have a very high velocity (> 700 km/s) 



reliably measured through an estimate of the distance obtained by ^ other methods than the 



dispersion measure. PSR B2224+65, po wering the Guitar Ne bula (iCordes et al. 



showing an X-ray tail simi 



km s ^ (IManchester et al. 



ar to Morla ( iHui fc Becker 



19931) and 



20071 ). has a speed of at least 865 



2005|). J0821-4300, the central compact o 



jject in the Puppis 



A supernova remnant has a speed possibly 



Winkler fc Petre 



2007; 



Becker et al. 



higher than 1000 km s ^ (IHui &: Becker 



2006 



20121 ). The spin-down power of PSR B2224+65 is 
similar to that of Morla {Erot=5.9 x 10'^^ erg s~^ and 1.2 x 10^'^ erg s~^), but the tails are 
different. While both tails have a similar brightness profile and show no spectral variations, 
the tail of PSR B2224+65 is misaligned (118°) with respect to the proper mo t ion o f the 



Bandiera 



(l2008h the 



pulsar. Alternative explanations have been proposed: according to 
feature could be produced by highly energetic particles which escaped the pulsar wind 
termination shock and are confined by an interstellar magnetic field. This scenario accounts 
for a faint counter-feature, detected on the opposite side of the pulsar w ith respect to the 



tail, b ut does not explain the absence of spectral variations along the tail (jjohnson fc Wang 



20101 ). Another possible s cenario is that the obse rved emission comes from a ballistic jet. 



similar to those of AGNs (Johnson fc Wang 



20101 ). but also this picture has some problems. 



The key to understand why these two pulsars do not emit through thermal bremsstrahlung 
lies in the ISM characteristics. This emission is detectable in X-rays only if both the ISM 
density and temperature are high enough (~0.2 atoms cm~'^ and ~10^-10^ K). It is hardly 
believable that the ISM around a CCO like J082 1-4300, lying inside its supernova remnant, 
could be denser than ~10~'^ atoms cm~^: this can explain the lack of a bremsstrahlung 
tail. In the case of PSR B2224+3205, the ISM density, estimated from the Ha bow-shock 
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emission to be < 1 cm~^ (IChatterjee fc Cordes II2002I ). is lower than the ISM density 
predicted around Morla. No constraints on the ISM temperature around PSR B2224+3205 
appear in the literature: it is not excluded that a low temperature of the ISM may 
contribute to the lack of a bremsstrahlung tail. 



Recently, 



Tomsick et al. 



(120 121 ) studied the case of the X-ray so urce IGR 



11014-6103, 



19681 ) ■ thus 



suggesting an association with the supernova remnant MSH 11-61 A (iKesteven 
inferring a velocity of 2400-2900 km s^^. IGR J11014-6103 shows a complex morphology, 
consisting of a point source, an elongated tail ~ 50" long, aligned to the direction of the 
proper motion, and a feature of ~4', nearly perpendicular to the proper motion. Spectral 
and spatial analyses suggest that IGR J11014-6103 could be a pulsar, therefore providing a 
natural explanation for the co-aligned feature as a bow-shock ram-pressure-confined pulsar 
wind nebula. Moreover, the feature resembles the tail of Morla: its luminosity peaks far 
from the pulsar and only a few photons are detected near the pulsar. The low statistics 
prevent us from studying further the similarities between the two tails. However, if the 
association with MSH 11-61A is confirmed, its distance would be 10 kpc, implying a tail of 
~ 2 pc, comparable to that of Morla. 



5.3. The Distance and inclination 



Our deep XMM-Newton observation of the field of J0357+3205 allowed us to better 



constrain the distance of the pulsar using the interstellar abso rption value. The G a. 



Dickey fc Lockman 



actic 



(11990f l 



absorption column in the direction of the pulsar, predicted by 
at ~500 pc based on the HI distribution, is fully compatible with our best fit value for 
the pulsar and the tail, pointing to a lower limit for the pulsar distance of ~300 pc. This 
result is in broad agreement with the value estimated by scaling the 7-ray fiux of the pulsar 



(jSaz Parkinson et al. 



2OIOI ): the method hinges upon the observed correlation between the 
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intrinsic 7-ray luminosity and spin-down power of the pulsar (lAbdo et al. Il2010a[). assuming 



a bea m correction factor of 1 for the 7-ray emission cone of all the pulsars (IWatters et al 



20091 ) ■ By applying this relation, we have a 7-ray efficiency of 0.33 and a distance of ~500 
pc. We obtain a new upper limit of ~900 pc by requiring the 7-ray efficiency to be less 
than 1. 

Another important piece of information could come from the serendipitous stars we found 
in the field: their distance estimates from optical analyses could provide lower limits on the 
distance of the pulsar. Future multi- filter optical observations will provide this limit. 
Under the hypothesis of a bremsstrahlung-emitting tail, we infer a pulsar velocity of ~1900 
km s~^. The estimated value of the pulsar proper motion, 165±60 mas yr~^, allows us to 
set another upper limit for the pulsar distance, being 2AtQ^]ipc for a purely transversal 
motion. However, as reported in section 15. 2[ a null inclination leads to unrealistic values 
(~10^ K) of the ISM temperature, so that an inclination angle larger than 60° is required. 
This lowers the upper limit to 1.21q3 kpc, still less constraining than that obtained by 
computing the pulsar efficiency. 

For a pulsar with a kick velocity of ~1900 km s^^, a proper motion of 165 mas yr~^ 
and a distance of 300 pc < -D < 900 pc, the inclination of the system with respect to the 
plane of the sky ranges between 68° and 83°, the density of the ISM between 0.15 and 0.3 
atoms cm~'^, and its temperature between 1 and 9 x 10^ K. The length of the tail is ~6.5 
pc, independently of the distance. 

Although the study of the kinematics of the pulsar is beyond the scope of this paper, 
the small component of the proper motion perpendicular to the Galactic plane points to 
a birth site outside the Galactic plane, suggesting that its pr ogenitor might have been a 



runaway star, as already introduced in iDe Luca et al. 



mm 
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6. Conclusions 

We have characterized the properties of the middle-aged 7-ray and X-ray pulsar 
J0357-I-3205 (Morla) and its tail. The X-ray spectrum of the pulsar is consistent with a 
magnetospheric non-thermal component plus thermal emission from hot spots; the lack of 
any detected thermal component from the whole surface makes Morla the coldest NS in its 
age range. 

Folding the X-ray photons at the expected pulsar period using the LAT ephemeris yields 
marginal evidence for modulation, with 2 peaks, neither of them aligned with the single 
gamma-ray peak (see Figure [3]). 

We confirmed the presence of a diffuse emission ~9' long, the largest tail of X-ray 
emission associated with any rotation-powered pulsar. Such an extended emission cannot 
be explained in terms of the usual bow-shock ram-pressure- dominated pulsar wind nebula. 
In fact, the following problems arise: 

- the existence of the tail is problematic for a pulsar with such a low E as Morla; 

- no bow shock has been resolved around the pulsar; 

- the lack of any nebular emission around the pulsar, where the particle acceleration is 
maximum, cannot be explained; 

- the nebular spectrum lacks any spatial evolution, a necessary signature of the radiative 
cooling of the electrons accelerated at the wind termination shock; 

- the very asymmetric brightness profile in the direction perpendicular to the main axis of 
the tail requires a large ad-hoc inhomogeneity of the ISM around the pulsar; 

- fitting the XMM-Newton nebular spectra with a powerlaw, we found the nebular Nh to be 
higher than the Galactic value and only marginally in agreement with the Nh of the pulsar. 

We propose a thermal bremsstrahlung model as an alternative explanation of the tail 
emission. In this scenario, the emission comes from the shocked ISM material heated up to 
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X-ray temperatures. This model gives full account of the peculiar features of the tail: 

- the lack of any detectable spatial evolution in the tail spectrum is due to the long 
bremsstrahlung cool-down time; 

- the peculiar asymmetries of the brightness profile can be interpreted in terms of small 
changes in the ISM density that strongly affect the brehmsstralung emissivity; 

- most of the energy in the pre-shock flow is carried by the ions, while the electron 
temperature is responsible for the X-ray emission; the Coulomb heating time of the 
electrons behind the shock is fully in agreement with the lack of any detected nebular 
emission near the pulsar; 

- the value of Nh measured in the tail agrees both with the value obtained for the pulsar 
and the Galactic value. 

This scenario allows us to estimate some parameters of the pulsar and of the ISM around 
it. For a bremsstrahlung-emitting tail we estimate the pulsar distance to be between 300 
and 900 pc. A pulsar velocity of ~1900 km s~^ is required, in agreement with the pulsar 
proper motion for distances of some hundreds parsecs and a high inclination. The mean 
density of the ISM is required to be ~0.2 atoms cm~^ and the temperature of some 10^ K. 
Th is type of hot gas us ually presents lower densities, but a denser phase (possibly detected 



by 



Yao fc Wang 



(120051 )) is predicted to be the result of the action of massive stellar winds 
and supernova explosions. However, the low metallicity we obtained from the spectral fit 
makes the explanation of this gas envelope more diffcult. 

A high inclination of the system (>70°) is predicted in both models. This value provides 
further support to the bremsstrahlung model, where the cool-down time demands for a 
longer tail, and adds to the problems of a synchrotron model. 

For all these reasons, we believe Morla's nebula to be the first example of a new 
"turtle's tail" class of thermally-emitting nebulae, ironically associated with fast moving 
pulsars. 
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Until now we have no clear evidence of other pulsars characterized by a tail emitting via 
thermal bremsstrahlung, possibly for the requirements of a very fast pulsar and a hot, 
dense ISM. Moreover, energetic pulsars can also generate classic synchrotron nebulae, that 
may outshine a bremsstrahlung component. The recently discovered nebula associated with 
the source IGR J11014-6103, suspected of being a high-velocity pulsar, resembles the tail 
of Morla: further studies of this feature could lead to a second example of a turtle's tail 
nebula. 

This work was supported by the contract Swift ASI-INAF 1-004-11-0. We warmly 
thank the ApJ referee for the useful comments and suggestions as well as the lASF-INAF 
staff who constructively participated in the discussion during my "Astrosiesta" presentation. 
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